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ABSTRACT 

Context. Detailed studies of the shape of dust emission spectra are possible thanks to the current instruments capable of simultaneous 
observations in several sub-millimetre bands (e.g., Herschel and Planck). The relationship between the observed spectra and the intrin- 
sic dust grain properties is known to be affected by the noise and the line-of-sight temperature variations. However, some controversy 
remains even on the basic effects resulting from the mixing of temperatures along the line-of-sight or within the instrument beam. 
Aims. Regarding the effect of temperature variations, previous studies have suggested either a positive or a negative correlation be- 
tween the colour temperature T c and the observed spectral index /fobs- Our aim is to show that both cases are possible and to determine 
the principal factors leading to either behaviour. 

Methods. We start by studying the behaviour of the sum of two or three modified black bodies at different temperatures. Then, with 
radiative transfer models of spherical clouds, we examine the probability distributions of the dust mass as a function of the physical 
dust temperature. With these results as a guideline, we examine the (T c ,/3 obs ) relations for different sets of clouds. 
Results. Even in the simple case of models consisting of two blackbodies at temperatures Tq and Tq + ATq, the correlation be- 
tween T c and fj obs can be either positive or negative. If one compares models where the temperature difference AT between the two 
blackbodies is varied, the correlation is negative. If the models differ in their mean temperature T rather than in AT , the correlation 
remains positive. Radiative transfer models show that externally heated clouds have different mean temperatures but the widths of 
their temperature distributions are rather similar. Thus, in observations of samples of such clouds the correlation between T c and /fobs 
is expected to be positive. The same result applies to clouds illuminated by external radiation fields of different intensity. For internally 
heated clouds a negative correlation is the more likely alternative. 

Conclusions. Previous studies of the (7c,/?) relation have been correct in that, depending on the cloud sample, both positive and neg- 
ative correlations are possible. For externally heated clouds the effect is opposite to the negative correlation seen in the observations. 
If the signal-to-noise ratio is high, the observed negative correlation could be explained by the temperature dependence of the dust 
optical properties but that intrinsic dependence could be even steeper than the observed one. 

Key words. ISM: clouds - Infrared: ISM - Radiative transfer - Submillimeter: ISM 



1. Introduction 

Thermal dust emission is increasingly important as a tracer of 
the dense interstellar clouds and of the star formation process. 
Earlier ground-based and satellite observations already showed 
that the long wavelength dust opacity is not constant and the 
variations are prob ably tracing grain coagulation and aggre- 
gation processes jC ambresv et alJ l200l1: del Burgo et al.1 t2003: 
Kram er et al. 2003; Lehtin en et alJ 12007). The sub-millimetre 
balloon borne experiments PRONAOS and Archeops enabled 
the study of the dust spectral index for a large number of in- 
terstellar clouds. By a coverage of the sub-millimetre regime, it 
was possible to start to separate the effects of the colour tem- 
perature, Tq, and of the dust emissivity spectral index, /?. The 
colour temperature Tc can be derived from the fit of a modified 
black body curve B(rc)v^ obs to the multi-wavelength observa- 
tions. The results indicated a negative correlation between the 
colour temperature and the spectral index with /fob s incre asing 
in the cold regions dDupac et al.ll2003b iDesert et al.1l2008l) . The 
results were not universally trusted because the variables are in- 
trinsically negatively correlated so that a small error in /Jobs can 
be compensated by a small error in Tc in the opposite direction. 
Therefore, and because of the other effects discussed in this pa- 
per, one must make the separation between the apparent spectral 
index /Jobs and the intrinsic spectral index (3 of the dust grains. In 



particular. lShettv et alJ ( f2009bl) and lShettv et al.1 (l2009al) showed 
that a negative correlation between Tc and /fobs could result 
from noise. However, the recent Planck and Herschel satellite 
results also show a decreasing /fobs(Tc) relation that, according 
to the authors, is too strong to be explained by the noise alone 
(lAnderson et al.l2010tlParadis et al.l2010tlVeneziani et al.l2010t 
iPlanck Collaboration et al.ll201 lblal) . The Herschel bands from 
100//m to 500 fim cover the peak of the dust emission. By cov- 
ering the longer wavelengths from 350 fim to millimetre waves, 
Planck is in principle in a better position to measure the emis- 
sion spectral index. However, for an accurate determination of 
/3 one also must constrain the temperature. For this reason the 
Planck data are being combined wit h far-infrared observations, 
for ex ample the 100/im IRAS data (IPlanck Collaboration et al.l 
1201 lbllah . 

The negative correlation of Tc and betao^ can be related 
to laboratory measurements of some interste llar dust analogues 
(lMennellaetal.lll998t ICoupeaud et al.1 120 111) and could be ex- 
plaine d by models developed for the emis sion of amorphous 
solids dMenv et alj|2007t iParadis et al.ll201 ll) . However, even if 
one accepts the negative correlation between /fobs and Tc as 
a real property of the observed radiation, there are still some 
hurdles before conclusions can be drawn regarding the intrin- 
sic properties of the dust grains. The mixing of different dust 
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temperatures along the line-of-sight, or more generally within 
the beam, means that the peak of the dust emission broadens 
and that the observed /fobs values are reduced. At the same time, 
the colour temperature will tend to overestimate the mass aver- 
aged physical dust temperature, possib ly resulting in a se rious 
underestimation of the cloud masses (Mal inen et al.|[20TTb . In 
this case the observed points would be displaced towards the 
lower right in the (7c, /fobs) plane. This has been interpreted 
as evidence that the negative correlation between the observed 
Tg and /fohs values a lso could r esult from temperatur e variations 
dShettv et al.ll2009al) . However. iMalinen et alTd201 lb carried out 
radiative transfer calculations for turbulent model clouds con- 
taining gravitationally bound cores. As long as the cores were 
externally heated, the result was an apparent positive correlation 
between Tc and /fobs- 

The purpose of this paper is to seek an explanation for this 
apparent contradiction and to establish what is the expected be- 
haviour for dense clouds and why. 

We start by describing our methods and the basic assump- 
tions in Sect. [2] The main results are presented in Sect. [3] We 
start with two and three layer models and, in this simple setting, 
examine the basic effects that the temperature mixing has on the 
observable Tc and /fobs values. In Sect. I3.3l we carry out radia- 
tive transfer modelling of a series of spherical model clouds to 
determine what kind of temperature probability distributions are 
expected for a collection of dense clouds. In Sect. l3.4l we use this 
information together with the two layer models to demonstrate 
the expected behaviour for quiescent clouds. Our final conclu- 
sions are presented in Sect. [4] 

2. The methods 

We examine the relations between the intrinsic dust temperature, 
T, and opacity spectral index, {3, with the corresponding param- 
eters derived from the analysis of the observed emission. These 
are the colour temperature, Tc, and the observed spectral index, 
/fobs, both deduced from the shape of the intensity spectrum. The 
observed parameters are calculated by fitting a modified black 
body fit, B v (7c)v 0Obs , to observations in the five Herschel bands 
at 100/im, 160//m, 250 //m, 350 yt/m, and 500 fim. We use the 
monochromatic intensity values and weighted least squares fits 
where the same relative uncertainty is assumed for all the bands. 

The simplest case where the observed spectral index /fobs can 
differ from the intrinsic f3 of the grains, is the two layer model. 
The source consists of two layers at different temperature both 
possibly still having the same intrinsic [3 value. We will also in- 
vestigate cases with three layers, i.e., mixtures of dust at three 
temperatures. The emission is assumed to be optically thin so 
that the observed intensities are a sum of the emission of the in- 
dividual layers. Thi s allows us t o comp are our findings directly 
with the results of IShettv et alj (|2009b). The model parameters 
that can be varied are the physical dust temperatures and the rel- 
ative column densities of the layers. The value of [3 is assumed 
to be 2.0 but we will also briefly look at cases with different f3 
values in the layers. 

We will ca rry out radiati ve transfer modelling of Bonnor- 
Ebert spheres (Bonnor 1956). These calculations provide first 
information on the temperature distributions that are likely to be 
found in real clouds. The density profiles are calculated for al- 
most critically stable configurations (stability parameter £=6.5) 
with a gas kinetic temperature of 10K. The precise values of 
these parameters are not crucial here because the dust temper- 
ature distributi ons depend mainly on the column density (e.g. 
Fischera 2011). We start radiative transfer modelling assuming 
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Fig. 1. The variation of the observed colour temperature and 
spectral index in two layer models (solid curves). The first layer 
is at fixed temperature of T = 10K. The other layer is at 1 1 K, 
12 K, or 14,K (black labels in the figure) and its relative mass 
is varied between 0.1 and 10.0 (four values are indicated on 
the curves). The intrinsic spectral index of both components is 
j8=2.0. 



that the clouds are heated by an external radiation field corre- 
spond ing to the interstellar radiation field (ISRF) model of Black 
d!994l) . The dust p roperties are those of the Milky Way dust 
with Ry-5.5 (Weingart ner & Draind lToOl). The dust tempera- 
ture profil es are derived with Mon t e Carlo radia tive transfer cal- 
culations dJuvela&Padoari 120031: Uuvelal 12005). Together with 
the density profiles, these can be converted to probability distri- 
butions of the column density as a function of dust temperature. 
The Bonnor-Ebert models are discussed in Sect. 13.31 

3. The results 

3.1. Two layer models 

Figure Q] shows the observed temperatures and spectral indices, 
Tc and /fobs, for a series of two layer models. The first dust layer 
has a temperature of T=10 K and the other one a temperature of 
1 1 K, 12 K, or 14 K. Both have an intrinsic opacity spectral index 
of {3=2.0. In Fig.Q~l the curves correspond to different mass ratios 
where the relative mass of the warmer dust is varied between 0. 1 
and 10. 

The correlation between Tc and /fobs is positive when the 
mass ratio of warm to cold dust is close to unity or slightly 
higher. When either temperature component is dominating the 
mass, the correlation between Tc and /fobs becomes negative. If 
one compares models with the same mass ratios but different 
temperature combinations (i.e., 10K+11K vs. 10K+14K) the 
correlation also is negative. 

A two layer model is too simple to be directly compared to 
actual observations. Nevertheless, one can appreciate the large 
range of possible (Tc, /fobs) combinations that are theoretically 
possible. In Fig.Q]each curve would end at a value of /fobs=2.0 
both when the mass ratio goes to zero or infinity. One notable 
feature is that for a large temperature difference, i.e., 10 K mixed 
with 14K, the curve passes through a region where Tc ex- 
ceeds the physical temperature of the both componen ts. Such 
behaviour was already observed by She ttv et aD ([2009b). 

In Figure [2] we use different intrinsic (3 values for 
the two layers with (3 =2.3, 2.2, 2.1, and 2.0, for the 
temperatures of 10 K, UK, 12 K, and 14 K, respectively. 
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Fig. 2. Modification of the two layer model of Fig.Q]with differ- 
ent intrinsic spectral indices. The 10 K cold component has an 
intrinsic spectral index of/? = 2.3. The warmer component is at 
a temperature of 11 K, 12 K, or 14 K (as indicated by the black 
labels) and has a spectral index j3 of 2.2, 2.1, and 2.0, respec- 
tively. The red circles correspond to these (T, f3) values of the 
individual dust components. 



This approximates the negative correlation reported in ob- 
servational studies (e.g. |5upac et al.l 120031: iDesert et al.| [2008; 
iPlanck Collaboration et all 1201 lblT Compared to Fig. [Q the 
curves are tilted to show stronger negative correlation between 
Tc and /fobs- The difference between /fobs and the average intrin- 
sic jS is smaller than in Fig. [TJ To some extent this is expected 
because the warm component, with the smaller f3, dominates 
the observed intensities. However, in Fig. |2]the combination of 
10K+14K has a minimum at /fobs ~1.55, less than 0.8 units be- 
low even the higher of the two /3 values. In Fig.[TJthe discrepancy 
between /fobs and the common f3 is larger. 



3.2. Three layer models 

We examine models with three temperature layers to make sure 
that the qualitative behaviour does not change as one takes the 
first step towards more continuous temperature distributions. 
Figure [3] shows the results for a mixture of dust at 8K, 10K, 
and 12 K. The starting point is a situation where the three tem- 
perature components have equal mass. The relative masses are 
then scaled between 0.1 and 1.0, one layer at a time. The first 
observation is that the /fobs values are now significantly lower 
because of the addition of the colder component at 8 K. 

When the column density of the 8K or of the 12 K compo- 
nent is varied, the behaviour is qualitatively similar to Fig. [TJ 
When the mass of the 10 K component is varied, the (Tc, /fobs) 
values fall on an almost straight line. The column densities were 
scaled in multiplicative steps of 1 .2 so that the distance between 
the plotted symbols corresponds to a constant step in log N. The 
corresponding displacements in the (Tc, /fobs) plane are quite 
regular, deviations becoming noticeable mainly when the rela- 
tive mass of the 8 K layer is modified. 

If one continues the curves for column density multipliers 
larger than 10, they will eventually reach the point with /fobs =2.0 
and Tc equal to the T of the dominant layer. The correlation 
between the Tc and /fobs values can again be either positive or 
negative, depending on which particular models are being com- 
pared. Thus the key question becomes, how the real clouds dif- 
fer from each other regarding the range of temperatures and the 
relative masses of the temperature components. To answer this 
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Fig. 3. As Fig. [TJbut for three layers at 8 K, 10 K and 12 K. The 
central point corresponds to identical mass in each temperature 
component and the curves are obtained by scaling the relative 
mass of one component at a time. For the component that is be- 
ing modified, the temperature is indicated in the figure (black 
labels) and the relative mass is shown for a few positions along 
the resulting curves (gray labels). All components have the same 
spectral index, /3 =2.0. 



question, we need models where the temperature distributions 
are solved self-consistently considering the balance of dust heat- 
ing and cooling. 



3.3. Bonnor-Ebert spheres 

We solve dust temperatures for a series of Bonnor-Ebert spheres 
with masses 0.2, 0.4, 0.8, and 1.6 Mq (see Sect. 13. The first 
frame of Fig. [4] shows the density and temperature profiles of 
the model clouds. The temperature at the surface of the clouds 
is 17-18 K and the central temperature decreases from ~12K 
in the 1.6 Mq cloud to ~8K in the 0.2 Mq cloud. The second 
frame shows the probability distributions of the column density 
on a line-of-sight through the cloud centre as a function of the 
temperature. The shape of the distributions is rather similar irre- 
spective of the cloud mass and the probability distributions are 
only shifted towards lower T as the column density of the model 
increases, i.e., when the mass of the Bonnor-Ebert sphere is de- 
creased. The situation is less clear when we consider the total 
dust mass. The main effect is still a shift in the temperature al- 
though the skewness of the distribution is larger for the clouds 
with the smallest mass. This suggests that the main variation for 
these externally heated clouds is in the mean temperature, not 
in the magnitude of the temperature variation nor in the relative 
mass of the (in relative terms) warm and cold components. 



3.4. The source of a positive correlation 

With the knowledge of the temperature distributions found in 
the radiative transfer models we can now return to the two layer 
models and use them to demonstrate the behaviour of an ensem- 
ble of clouds. The models consist of two layers at temperatures 
To and To + AT. Because we want to compare these results to 
the spectra calculated for the Bonnor-Ebert models and want to 
avoid the complication of the contribution of the stochastically 
heated dust grains to the 100/im intensity, these fits were done 
without the 100 //m points. Figure [3] shows the results for differ- 
ent values of AT, for a set of models with 7q=8.0, 10.0, 12.0, 
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Fig. 4. Frame a: The radial density profiles (solid curves and the 
scale on the left) and the temperature profiles (dashed curves and 
the scale on the right) for Bonnor-Ebert spheres of 0.2, 0.4, 0.8, 
and 1.6 Mq (curves from left to right). Frame b: The probability 
distributions of the column density as the function of the dust 
temperature for a line-of-sight through the centre of the cloud 
(solid lines). The curves correspond to the 0.2, 0.4, 0.8, and 
1 .6 Mq models, from left to right. The dashed lines show the 
corresponding probability distributions for the total dust mass 
of the model cloud. All the probability distributions have been 
normalised to a maximum value of 1.0. 



or 14.0 K. A ratio of 4:1 is assumed for the column densities of 
the cold and warm components. The same figure also shows the 
results for the Bonnor-Ebert spheres. In the two cases, the clouds 
are illuminated either by the full ISRF or by the ISRF attenuated 
by an external layer of dust with Ay = 2 m . 

When AT is kept constant but Tq is varied, the behaviours 
are qualitatively similar to the Bonnor-Ebert spheres in that the 
correlation between Tq and /fobs is positive. On the other hand, if 
Ar is varied keeping Tq constant, the correlation becomes nega- 
tive. The calculations were repeated with column density ratios 
4:1, 3:1, 2:1, and 1:1 for the T Q values of 8, 10, 12, and 14K, 
respectively (see Fig. [5] the open symbols). This does not have 
a major impact on the correlations between the Tc and the /fobs 
parameters. 

This key conclusion is that the sign of the correlation entirely 
depends on the models o ne chooses to comp are. Therefore, it is 
incorrect to interpret the IShettv et~ai] (l2009ah results as a gen- 
eral proof that the temperature variations always cause a nega- 
tive correlation between Tq and /fobs- 
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4. Discussion 



In their article She ttv et alj d2009bl) examined the effect of noise 
and of line-of-sight temperature variations on the observed T 



Fig. 5. The dependence between Tc and /? for a model of two 
layers at temperatures To and Tq + AT, with a ratio 4: 1 between 
the column densities. Each solid curve corresponds to a single 
value of Ar, with 7o=8, 10, 12, and 14 K along the dashed lines. 
The open symbols show corresponding results for column den- 
sity ratios 4:1, 3:1, 2:1, and 1:1 for the T values of 8, 10, 12, 
and 14 K, respectively. The two dotted lines and the square sym- 
bols show the values for Bonnor-Ebert spheres of 0.2, 0.4, 0.8, 
and 1.60 Mq) with no external shielding or spheres shielding by 
Ay = 2 m . In this figure, the 100 /im intensities were not used. 



and ft values. They concluded that both these factors can pro- 
duce a negative correlation between /fobs and Tq or at least sig- 
nificantly affect t he observed /fo b s 7c re lations. Contrary results 
were obtained by Malinen et al] (1201 ll) who analysed synthetic 
sub-millimetre observations. The cloud model was based on 
magnetohydrodynamic simulations where the self-gravity had 
produced a number of cores, some with very high column den- 
sities. Radiative transfer modelling was carried out and the syn- 
thetic observations were analysed. The /fobs values were always 
below the intrinsic (3 but a clear positive correlation was seen 
at low temperatures, i.e., associated with the set of dense cores. 
The proposed explanation was that one is again mixing cold and 
warm dust. In Malinen et al.l(l201ll) the correlations were derived 
for all pixels in the synthetic maps. We present in Appendix lAl 
further analysis of those data, concentrating on the locations of 
the dense cores. 

An explanation based on the temperature mixing is, however, 
incomplete because we have seen that the result depends on the 
way how the warm and the cold components are mixed. Most 
importantly, one needs to know what is the basic difference of 
the clouds one is observing. The radiative transfer calculations 
carried out with spherical model clouds (Sect. 13 . 3b showed that, 
in the first approximation, the main difference is their mean tem- 
perature. The shape of the probability distribution for the col- 
umn density as a function of temperature remained rather similar 
and the distribution was only shifted along the temperature axis 
(Fig. |4j5). In this scenario, the simplest description of an ensem- 
ble of clouds is a set of two layer models with temperatures Tq 
and Tq + AT where the clouds differ in Tq but not in Ar. For un- 
resolved clouds the changes in the shape of the N(T) distribution 
were more noticeable but still not sufficient to alter the conclu- 
sion. The expected behaviour is a positive correlation between 
T c and /fobs (Fig. 01. 

If the parameter Ar is varied instead, the result is a strong 
nega tive correlation betw een the temperature and the spectral in- 
dex. HhettxetaD (]2009b) did not directly show the T c (J3 0bs ) re- 
lations as a function of model parameters. Comparing their Figs. 
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Fig. 6. The parameters Tc and f3 for a set of internally heated 
Bonnor-Ebert models with masses of 0.4, 0.8, 1.6, 3.2, and 
6.4 Mq. The solid lines connect models with an internal heat- 
ing source of given luminosity, 0.0, 0.1, 0.3, 1.0, or 3.0 solar lu- 
minosities. The dashed lines connect models of the same mass. 
The parameters have been derived assuming the source in unre- 
solved. 



3 and 4 one can find either positive or negative correlation. For 
example, the correlation is negative between the 10K+15 K and 
the 10K+20K cases, precisely when the parameter AT is be- 
ing changed. In Malinen et al] (1201 ll) . the correlation between 
Tc and /fobs also became predominantly negative when internal 
heating was applied to the sources. This is schematically consis- 
tent with the above picture. When cores are heated by sources of 
different luminosity, the temperature gradients will be different. 
In other words, one is primarily modifying the Ar parameter. In 
Fig- El this would correspond to the cores being distributed in a 
direc tion perpendicular t o the lines of constant Ar (see also Fig. 
8 of lShettvetafll200 9b]). To illustrate this further, Fig. shows 
results for a set of Bonnor-Ebert models with internal heating. 
The heating source is a 5800 K black body that is placed at the 
centre of the Bonnor-Ebert sphere. The figure shows the results 
for the cores assuming they are not resolved in the observations. 
As the luminosity of the point source increases, the correlation 
between Tc and /fobs becomes negative. This is true whether one 
is comparing models of different mass or models with different 
source luminosity. 

For large surveys (e.g. iPlanck Co llaboration et all 1201 lbh . 
one further consideration is the different intensity of the radi- 
ation field that is heating each cloud. Figure [5] already showed 
two cases, one for the normal ISRF, and one for the ISRF ex- 
tincted by Ay = 2 m . To extend this test to higher radiation field 
intensities, Fig. [7] shows the (7c,/fobs) values for the Bonnor- 
Ebert models when the radiation field intensity is multiplied by 
factors ^=1.0, 2.0, 4.0, and 8.0. The correlation between Tc 
and /fobs is positive whether one is comparing clouds of dif- 
ferent mass or clouds of given mass but subjected to different 
levels of isotropic external radiation field. If the external heating 
is not isotropic, the situation corresponds qualitatively to two 
layer models where the temperature difference between the lay- 
ers, Ar, is changing from source to source. This could again lead 
in observations to a negative correlation between the temperature 
and the spectral index. 

The result of the above discussion is that if the real f3(T) rela- 
tion is flat, the observed /fob s (7c) relations will be very different 
depending on the set of sources examined. For clouds heated 
mainly by an external isotropic radiation field, the observed /fobs 
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Fig. 7. (Tc, /fobs) values in synthetic observations of Bonnor- 
Ebert spheres illuminated by external radiation fields of different 
intensity. The cloud masses are 0.2, 0.4, 0.8, and 1.6 Mq and the 
ISRF is scaled by factors x= 1 -0, 2.0, 4.0, and 8.0, as indicated 
in the figure. The open circles show the case when the ISRF 
0tf=1.0) is attenuated by an external layer with Ay = 2.0™ (same 
as in Fig. [5]l. 

should be an increasing function of the colour temperature, i.e., 
the correlation between Tc and /fobs is positive. When internal 
or anisotropic external heating dominates, as in very active star 
forming regions, /fob s (7c) should be a decreasing function of 
temperature, i.e., the correlation is negative. These are merely 
the qualitative conclusions and detailed modelling is needed to 
estimate how the f3(T) relation may be modified by the temper- 
ature mixing. Only when models are constructed for individual 
sources does it become possible to take into account all the rele- 
vant factors like the precise source geometry, the anisotropies of 
the local radiation field, and the location and luminosity of the 
relevant point sources. 

5. Conclusions 

Even when the grain optical properties are independent of the 
temperature, the line-of-sight temperature variations can cause a 
correlation between the colour temperature Tc and the observed 
spectral index /fobs that can be either positive or negative. The 
sign of the correlation depends on the nature of the clouds whose 
physical differences lead to the dispersion of the observed Tc 
and /fobs values. Exam ination of a set of sph erical model clouds 
confirmed the result of Mali nen et al.l {201 1), a positive correla- 
tion between Tc and /fobs in the case of externally heated clouds. 
The correlation remains positive when the clouds are heated by 
isotropic radiation fields of different intensity. However, internal 
heating sources and anisotropic external heating can make the 
correlation negative. Quantitative estimates cannot be derived 
without detailed modelling and the knowledge of the exact na- 
ture of the sources observed. 
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Appendix A: Analysis of clumps in a MHD model 
cloud 



Mali nen et alj d201 ll) presented the analysis of synthetic obser- 
vations of a model cloud obtained from a 3D MHD run with 
self-gravity. The model cloud had a linear size of lOpc and the 
resolution of the maps was ~0.005pc. Their Fig. 15 shows the 
(7c, /fobs) relations for all pixels in the synthetic maps. The re- 
lation is shown both when the model is heated only by an exter- 
nal radiat ion that corres ponds to the local interstellar radiation 
field (Ma this et al.ll9 83) and when internal heating sources were 
added to gravitationally bound cores. The masses of the cores are 
in the range of ~ 1-39Mq. The luminosities of the added internal 
sources were determined by the core mass and are in the range 
of ~0.3-12 solar lumino s ities. For the details of those calcu- 
lations, see iMalinen et alj d201 lb . The colour temperatures and 
spectral indices were based on modified black body fits to data 
at 100/im, 160/itn, 250 fim, 350 fim, and 500 fim, with the data 
convolved with a gaussian with FWHM equal to 20 pixels. 

We repeat the analysis using the same synthetic surface 
brightness observations but concentrating on the locations of 
the gravitationally bound cores. The model contains 40 such 
cores but internal radiation sources were added only in 34 cores 
that were sufficiently resolved in the automatic mesh refinement 
(AMR) calculations. Therefore, we concentrate on the (7c, /fobs) 
values towards those 34 cores, before and after the addition of 
the interna l heating. 

Unlike Malinen et al 1 (1201 lh . we examine the data for three 
orthogonal viewing direction. Also, we only use wavelengths 
from 160//m to 500 //m. The lOO/vm is affected by the small 
grain emission and is often omitted for the analysis of the large 
grains. We assume a cloud distance of 500 pc or 1 .0 kpc and anal- 
yse the surface brightness values obtained for a single beam to- 
wards each of the cores. The posit ions are known fro m the anal- 
ysis of the 3D density cube (see Malinen et al. 201 1). The beam 
size is 40" that, depending on the distance, corresponds to 20 or 
40 pixels. 

The derived spectral index and colour temperature values 
are shown in Fig. I A. 1 1 These are in qualitative agreement with 
the results we obtained for spherical models. For externally 
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Fig.A.l. The observed spectral indices and colour tempera- 
tures towards t he se lf-gravitating cores in a model presented in 
Mali nen et al.1 d201 ll) . The frames a and b show the situation 
without and the frames c and d with heating sources inside the 
cores. The assumed cloud distance is 500 pc in frames a and c 
and 1 kpc in the frames b and d. The symbol sizes correspond to 
the mass of the core (~1-39Mq). The surface brightness maps 
were calculated for three orthogonal directions. For each core, 
the lines (of random colour) connect the values obtained for 
three different directions of observation. The linear correlation 
coefficient between Tc and /fobs (including all the plotted points) 
is given in the lower right hand corner of each frame. 



heated cores (frames a and b) the correlation between Tq and 
/fobs is positive. The correla tion is very weak but consistent 
with the Mali nen et alj d201 lh plots that contained data for the 
whole extent of the model cloud. The correlation coefficients 
are shown in the figure. With internal heating sources (frames 
c and d), the correlation between Tq and /fobs is clearly negative. 
Furthermore, the relation appears to be no n-linear and reminis - 
cent of the error bananas produced by noise Shett y etalJ d2009a). 
However, no noise was added to these surface brightness values 
and the scatter in Fig. IA. H is caused entirely by the real temper- 
ature variations. The (Tc, /fobs) values do not show very signif- 
icant differences between the viewing directions (i.e., the trian- 
gles in Fig. I A. 1 1 are relatively small). This is natural when the 
emission is optically thin and most of the emission originates in 
a single core along the line of sight. If some lines of sight had 
crossed cores of different temperature, more significant variation 
would have been observed. Similarly, at least for this particular 
model, the results are not sensitive to the resolution of the obser- 
vations. 
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